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Abstract 

Background: Fungi are asexually and sexually reproducing organisms that can combine the evolutionary 
advantages of the two reproductive modes. However, for many fungi the sexual cycle has never been observed in 
the field or in vitro and it remains unclear whether sexual reproduction is absent or cryptic. Nevertheless, there are 
indirect approaches to assess the occurrence of sex in a species, such as population studies, expression analysis of 
genes involved in mating processes and analysis of their selective constraints. The members of the Phiolocephala 
fortinii s. I. - Acephala applonoto species complex (PAC) are ascomycetes and the predominant dark septate 
endophytes that colonize woody plant roots. Despite their abundance in many ecosystems of the northern 
hemisphere, no sexual state has been identified to date and little is known about their reproductive biology, and 
how it shaped their evolutionary history and contributes to their ecological role in forest ecosystems. We therefore 
aimed at assessing the importance of sexual reproduction by indirect approaches that included molecular analyses 
of the mating type {MAT) genes involved in reproductive processes. 

Results: The study included 19 PAC species and > 3, 000 strains that represented populations from different hosts, 
continents and ecosystems. Whereas A opplonoto had a homothallic (self-fertile) MAT locus structure, all other 
species were structurally heterothallic (self-sterile). Compatible mating types were observed to co-occur more 
frequently than expected by chance. Moreover, in > 80% of the populations a 1:1 mating type ratio and gametic 
equilibrium were found. MAT genes were shown to evolve under strong purifying selection. 

Conclusions: The signature of sex was found in worldwide populations of PAC species and functionality of MAT 
genes is likely preserved by purifying selection. We hypothesize that cryptic sex regularely occurs in the PAC and 
that further field studies and in vitro crosses will lead to the discovery of the sexual state. Although structurally 
heterothallic species prevail, it cannot be excluded that homothallism represents the ancestral breeding system in 
the PAC. 



Background 

The origin and maintenance of sexual reproduction is a 
controversially discussed topic in evolutionary biology as 
reflected by the multitude of theories that have been 
proposed to explain why sexual reproduction, although 
highly costly, is widely occurring in nature [1-6]. Asex- 
ual organisms have a two-fold advantage over sexual 
conspecifics and can effectively disseminate [7-10]. In 
contrast, sexual reproduction efficiently eliminates dele- 
terious mutations [11] and creates genetic variation that 
favors natural selection and accelerates adaptation to 
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changing environments [1]. However, many species are 
capable of reproducing both sexually and asexually and 
illustrate how difficult it is to provide a general explana- 
tion on the evolutionary significance of sex. Fungi com- 
bine the advantages of the two reproductive modes. 
Several reasons were proposed for why the cost of sex 
compared to asex is lower in fungi than for animals and 
plants [12] because (i) fungi can be isogamous and thus 
the contribution of resources to the zygote by the 
gametes is limited and (ii) many fungi are also homo- 
thallic (self-fertile) and do not depend on finding a com- 
patible mate which reduces the cost of sex, whereas 
others are heterothallic (self-sterile) and mating is regu- 
lated by mating type factors, (iii) Moreover, the majority 
of fungi can alternate between asexual and sexual 
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reproduction and thus sexuality can be adjusted to when 
opportunity costs are low, for example at the end of the 
growing season of a host plant on which a fungus is 
dependent when adverse conditions are disadvantageous 
for somatic growth. Sex is also linked to essential pro- 
cesses such as the formation of resistant spores that are 
able to survive unfavourable conditions and enable new 
genotypes to be spread into new environments. Once 
the link between sex and such essential processes has 
evolved, selection against frequent sexual recombination 
might be less effective. 

However, for many fungi, especially for filamentous 
ascomycetes, only part of their life cycle is known. 
These taxa are classified as Deuteromycota or "fungi 
imperfecti" due to the lack of sexual morphology [13], 
but it is unclear whether sexual reproduction is absent, 
rare or cryptic because sexual morphology is often diffi- 
cult to observe in nature or in the laboratory [14,15]. 
Thus, the importance of sexual reproduction in natural 
populations of such species remains an open question 
which can be addressed by direct and indirect 
approaches [16]. 

The direct apporach consists in searching for the sex- 
ual state (teleomorph) in the field or in the laboratory. 
However, it is often difficult to induce the teleomorph 
in vitro as many factors (e.g. nutrient media, tempera- 
ture, light exposure, selection of compatible mating 
types) need to be optimized for a successful induction 
[17,18]. Moreover, although sexual structures obtained 
in the laboratory indicate that the ability for sexual 
reproduction has not been lost, its importance in the 
field remains to be established by population studies 
and by monitoring the teleomorph in the field [16]. 
Indirect approaches comprise population studies that 
test the null hypothesis of random mating and include 
analysis of gametic disequilibrium, mating type ratios, 
genotypic diversity and phylogenetic analysis [19,20]. 
Moreover, functionality of genes involved in mating pro- 
cesses can be tested by expression analyses and provide 
further evidence for sexual reproduction [21]. Alterna- 
tively, estimates of selective pressures acting on such 
genes might indicate whether functionality is preserved 
[22-26]. 

In this study we aimed at elucidating the importance 
of sexual reproduction in the Phialocephala fortinii 
sensu lato - Acephala applanata species complex (PAC) 
that belongs to the dark septate endophytes, a polyphy- 
letic group of ascomycetes with characteristic melanized, 
septate hyphae that commonly colonize roots of woody 
plant species [27]. The PAC is composed of more than 
20 species, eigth of which were formaly described 
[28,29], and they dominate the endophytic assemblages 
in roots of conifers and members of the Ericaceae in the 
northern hemisphere from polar to subtropical regions 



[30-32]. PAC species form communities of up to ten 
species [32,33] and species abundance distributions 
within these communities ususally follow a hyperbolic 
distribution with a few abundant species and many 
"rare" species, consistent with the community structures 
of many other organismal groups [34]. Interestingly, no 
biogeographical pattern was found for PAC species and 
community compositions did neither correlate with 
composition of the host species nor with climate [32] 
supporting the the hypothesis that "everything is every- 
where" [35]. Despite the broad geographical occurrence 
of this fungal species complex, the teleomorph has 
never been observed for any of its species. 

In ascomycetes, the mating type (MAT) locus exercises 
key regulatory functions involved in mating processes 
and usually defines homothallism and heterothallism 
[36,37]. In heterothallic ascomycetes the MAT locus is 
characterized by two alternative forms, MAT 1-1 and 
MAT 7-2, called idiomorphs because they contain dis- 
similar DNA sequences at the same chromosomal loca- 
tion. These encode proteins with conserved DNA 
binding domains that are involved in transcription regu- 
lation leading to the attraction of compatible mating 
types in order to initiate the mating process [37] and 
are also involved in controlling the regulation of inter- 
nuclear recognition in later steps of sexual development 
[38-40]. In contrast, in homothallic species, a single 
individual generally contains all of the MAT genes and 
is thus capable of selfing [41]. In filamentous ascomy- 
cetes the simpliest idiomorphs consist of MAT 1-1 that 
carries one gene called MAT1-1-1 encoding a protein 
with an a -domain binding motif and MAT 1-2 that con- 
tains one gene called MAT 1-2-1 encoding a protein 
with a high mobility group (HMG)-binding domain 
[42,43]. However, additional MA T genes are regularely 
found in ascomycetous species [42]. For example, in the 
helotialean species Pyrenopeziza brassicae and Rhynchos- 
porium secalis that are closely related to the PAC [44] 
up to two additional MA T genes have been identified in 
the MAT1-1 idiomorph. One is designated MAT1-1-3, 
encodes a protein with an HMG-binding domain and is 
present in both species, whereas MAT 1-1-4 encodes a 
metallothionein-like protein and is only present in P. 
brassicae [45,46]. 

In a recent study we cloned the MAT genes of eight 
PAC species [44] . and expected the MA T idiomorphs to 
be similar to those of their helotialean relatives. Seven 
examined species showed a heterothallic organization of 
the MAT locus because strains either contained the 
MAT 1-1 idiomorph carrying MAT 1-1-1 and MAT 1-1-3 
or the MAT 1-2 idiomorph carrying MAT 1-2-1. In con- 
trast, in A. applanata the MAT locus structure was indi- 
cative of homothallism because all of the three MAT 
genes were present in single strains. However, it was 
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neither possible to induce sexual reproductive structures 
in crossing experiments nor to show that MAT genes 
were expressed under laboratory conditions [44]. In con- 
trast, analysis of multi-locus gametic disequilibrium for 
a limited number of Swiss populations of four PAC spe- 
cies showed that in most populations the index of asso- 
ciation I A [47] did not deviate significantly from zero 
which is indicative of recombination [33,48]. Because 
PAC species are successful root colonizers of a multi- 
tude of plant species in many ecosystems, we hypotesize 
that sexual reproduction occurrs as it confers the ability 
to successfully adapt to different environments 
worldwide. 

In the present study we aimed at gaining a more com- 
plete picture about the importance of sexual reproduc- 
tion in this species complex. In particular, we determined 
(i) the MAT locus structure for 11 additional PAC spe- 
cies, (ii) whether opposite mating types in populations of 
the PAC species deviated from the 1:1 ratio expected 
under random mating, (iii) the deviation from gametic 
equilibrium in the fungal populations, (iv) the spatial dis- 
tribution of mating types at the collection sites, and (v) 
the selective pressures acting on MAT genes. 

Results 

Organization of the MAT locus in PAC species 

The multiplex PCR amplified either a MAT 1-1 specific 
fragment of -550 bp or a MAT1-2 specific fragment of 
-750 bp in strains of any PAC species, except in the 
homothallic A. applanata for which no fragments were 
amplified, as expected due to its MAT locus organiza- 
tion [44] (Figure 1). During the screening of the collec- 
tions, occasionally (in < 2.5% of the screened strains) 
both fragments were amplified in single strains, indicat- 
ing potential homothallism, and in < 1.0% of the strains 
an amplification failure of the PCR product was 
recorded. However, re-analyzing newly prepared single- 
hyphal-tip cultures from a subset of eight of these 
strains, an amplification of single idiomorph-specific 
fragments was always obtained (see Additional file 1). 



The MAT locus of 11 additional PAC species was 
sequenced and characterized for one MAT 1-1 and 
MAT 1-2 strain per species (Table 1), and its structure 
was found to be congruent with the MAT locus struc- 
ture of the seven previously studied species [44]. The 
MAT 1-1 idiomorphs contained the MAT1-1-1 and 
MAT 1-1-3 gene whereas the MAT 1-2 idiomorphs 
included the MAT 1-2-1 gene. 

Mating type ratios and gametic disequilibrium in PAC 
populations 

The vast majority of populations (> 80%) did not deviate 
significantly from the expected 1:1 ratio in mating types 
(Table 2). After applying Bonferroni correction to adjust 
for errors of multiple comparisons (type 1 error), no sig- 
nificant deviations remained. Moreover, even those 
populations with < 10 individuals contained both mating 
types. Only a single population was found that did not 
include both mating types. Similarly, only six out of 52 
populations showed significant gametic disequilibrium 
and these were restricted to populations of P. subalpina 
(Table 2), Furthermore, in only two of these P. subal- 
pina populations both unequal mating type frequencies 
and gametic disequilibria were recorded. 

Spatial distribution of mating types in selected study 
sites 

The spatial distribution of mating types was analyzed for 
four PAC species in ten Swiss populations. In each 
population, strains of opposite mating type were found 
distributed over the whole study site (see Additional file 
2). Strains of different mating types were regularly iso- 
lated from the same grid points indicating that they can 
be found in close physical proximity (Table 3). In addi- 
tion, measures of association [Q] between opposite mat- 
ing types for significant spatial structures ranged from 
weakly negative to strongly positive (Table 3), meaning 
that co-occurrence of strains with different mating types 
is often more frequently observed than expected by 
chance. 




1618 bp 
1018 bp 

500 bp 

Figure 1 Example of multiplex PCR. Multiplex PCR amplification of idiomorph specific bands in selected strains of different PAC species, using 
the primers Pf_HMG_R.03, Pf_HMG_F4 and Pf_MAT1-1F1c. 
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Table 1 Strains for which the complete mating type idiomorph was sequenced 


Species 


Strain 


Idiomorph 


Geographic origin 


nost/oUDSiraie 


Reference/ 
Collector 


GenBank accession 
No. 


Phialocephola 
turicensis 


1_120_2 


MAT1-1 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347273] 




1_1 17_3 


MAT1-2 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347292] 


Phialocephola letzii 


2_120_3 


MATI-I 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347274] 




2_1 16_2 


MAT! -2 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347293] 


Phialocephola 
europoeo 


3_136_1 


MATI-I 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347275] 




3_122_3 


MAT! -2 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347294] 


Phialocepholo 
helvetica 


4_140_4 


MATI-I 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347276] 




4_153_2 


MAT1-2 


Zurichberg, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347295] 


Phialocepholo 
uotolensis 


5_265_3 


MATI-I 


Uetliberg; Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347277] 




5_265_4 


MAT1-2 


Uetliberg; Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347296] 


Phialocephola 
subalpino 


6 30 4 


MAT1-1 


Bodmeren, Switzerland 


root of Picea abies 


Zaffarano et al. 
(2010) 


[GenBank:HM347278] 




6 37 6v 


MAT1-2 


Bodmeren Switzerland 


root of Vaccinium 
myrtillus 


Zaffarano et al. 
(2010) 


[GenBank:HM347297] 


Phialocephola fortinii 
s.s. 


7 K93 444 a 


MAT1-1 


Suonenjoki, Finland 


rnrvK nf Pinu<\ svlvptfris 


z_a 1 1 a i ai i \j a i . 

(2010) 


[GenBank:HM347279] 




7_6_7v 


MAT! -2 


Bodmeren, Switzerland 


root of Vaccinium 
myrtillus 


Zaffarano et al. 
(2010) 


[GenBank:HM347298] 


CSP8 


8_L_05-7 


MATI-I 


Sheep Creek, USA 


n.a. 


J.B. Zaerr 


[GenBank:HM347280] 




8_L_02-5 


MATI-2 


Sheep Creek, USA 


n.a. 


J.B. Zaerr 


[GenBank:HM347299] 


CSP9 


9_ST_43-1 


MATI-I 


Walder Steinenberg, 
Switzerland 


root of Picea abies 


N. Brenn 


[GenBank:HM347281] 




9_ST_02-3 


MATI-2 


Walder Steinenberg 
Switzerland 


root of Picea abies 


N. Brenn 


[GenBank:HM347300] 


CSP10 


10_CS_15L_2 


MATI-I 


Court, Switzerland 


root of Picea abies 


C.R. Grunig 


[GenBank:HM347282] 




1 0_CS_9_1 


MATI-2 


Court, Switzerland 


root of Picea abies 


C.R. Grunig 


[GenBank:HM347301] 


CSP11 


11_AF196-16 


MATI-I 


Varena plantation, Lithuania 


root of Picea abies 


A. Menkis 


[GenBank:HM347283] 




11_AF200-1 


MATI-2 


Varena plantation, Lithuania 


root of Picea abies 


A. Menkis 


[GenBank:HM347302] 


CSP12 


12_405_1 


MATI-I 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347284] 




12_418_1 


MATI-2 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347303] 


CSP13 


13_422_6v 


MATI-I 


Creux du Van, Switzerland 


root of Vaccinium vitis- 
ideae 


V. Queloz 


[GenBank:HM347285] 




13_414_2 


MATI-2 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347304] 


CSP14 


14_426_4 


MATI-I 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347286] 




14_444_6v 


MATI-2 


Creux du Van, Switzerland 


root of Vaccinium 
myrtillus 


V. Queloz 


[GenBank:HM347305] 


Ljr 1 J 


1 D A UZ J 


IVIri III 


iNUUllUdy, Ujn 


n.a. 


I Mill 
J. I 1 1 1 1 


r C p n R a n V ■ \-\ h A 3 A 7 1 9, 7 1 
[ocl IDdl IK.nivl 3^x1 Zo/ J 




15_A_08-4 


MATI-2 


Noonday, USA 


n.a. 


J. Hill 


[GenBank:HM347306] 


CSP16 


16_428_6v 


MATI-I 


Creux du Van, Switzerland 


root of Vaccinium vitis- 
ideae 


V. Queloz 


[GenBank:HM347288] 




16_404_2 


MATI-2 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347307] 


CSP18 


18_410_1 


MATI-I 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347290] 




18_439_1 


MATI-2 


Creux du Van, Switzerland 


root of Picea abies 


V. Queloz 


[GenBank:HM347308] 


CSP19 


19_1 559-1 


MATI-I 


Parco del Pollino, Italy 


root of Abies alba 


V. Queloz 


[GenBank:HM347291] 




19_1 545-1 


MATI-2 


Parco del Pollino, Italy 


root of Abies alba 


V. Queloz 


[GenBank:HM347309] 



1 n.a. = information not available 
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Table 2 Mating type ratios for populations of species belonging to the Phialocephala fortinii s.l. - Acephala applanata 
species complex (PAC) 



Species Study site 



Country 



Marker 
type 1 



Number of 
strains 



Mating type ratios 



Gametic 
disequilibrium 













Number of 
MLH 


MAT1- 
1 


MAT1- 
2 


Significance 2 Number of 
MLH 3 


I 4 
'a 


P. 


Rothwald 


Austria 


M 


24 


7 


4 


3 


n.a. 




turicensis 






















Etang de la 


Switzerland 


M 


203 


12 


9 


3 


* 12 


0.28 




Gruere 




















Uetliberg 


Switzerland 


RFLP 


145 


17 


9 


8 


17 


0.07 




Zurichberg 


Switzerland 


RFLP 


48 


13 


7 


6 


13 


0.13 


P. letzii 


Parabock 


Switzerland 


M 


58 


1 1 


6 


5 


1 1 


0.69 




Rothwald 


Austria 


M 


46 


16 


7 


9 


16 


-0.01 




Zurichberg 


Switzerland 


RFLP 


52 


18 


13 


5 


* 18 


0.12 


P. 


Redwater 


Canada 


M 


19 


7 


1 


6 


n.a. 




europoea 






















Alaska 


USA 


M 


19 


12 


7 


5 


13 


0.41 




Catherine Creek 


USA 


M 


17 


8 


5 


3 


n.a. 






Sierra Nevada 


USA 


M 


10 


9 


7 


2 


n.a. 






Bodmeren 


Switzerland 


M 


28 


19 


7 


12 


17 


-0.18 




Scatle 


Switzerland 


M 


36 


31 


12 


19 


25 


0.02 




Rothwald 


Austria 


M 


86 


40 


23 


17 


38 


-0.12 




Alpes Maritimes 


France 


M 


77 


43 


24 


19 


34 


0.05 




Creux du Van 


Switzerland 


M 


123 


51 


27 


24 


43 


-0.08 




Ucl Ucl cl ILc 


OVVI LZcl Idl IU 


IVI 


1 7Q 
I / y 




19. 
zo 


3^ 

jO 


A9 


-n n^ 




Runcaglia 


Switzerland 


M 


95 


45 


22 


23 


41 


0.02 




Zurichberg 


Switzerland 


RFLP 


130 


46 


26 


20 


46 


0.04 


P. 


Yuba 


USA 


M 


25 


5 


3 


2 


n.a. 




helvetica 






















Ala Archa 


Kirgistan 


M 


18 


7 


7 


0 


n.a. 






Baikalsee 


Russia 


M 


29 


11 


4 


7 


12 


0.16 




Mansfield 


Canada 


M 


63 


16 


4 


12 


16 


0.46 




Sierra Nevada 


USA 


M 


47 


30 


17 


13 


30 


0 




Zurichberg 


Switzerland 


RFLP 


40 


15 


11 


4 


15 


-0.03 




Runcaglia 


Switzerland 


M 


98 


36 


20 


16 


36 


0.02 


P. 


Alpes Maritimes 


France 


M 


17 


6 


4 


2 


n.a. 




uotolensis 




















P. 


Creux du Van 


Switzerland 


M 


13 


8 


6 


2 


n.a. 




subolpina 






















Etang de la 


Switzerland 


M 


1 5 


14 


8 


6 


14 


0.05 




Gruere 




















Massif Central/La 


France 


M 


70 


30 


14 


16 


29 


-0.08 




Fage 




















Sheep Creek 


USA 


M 


12 


7 


2 


5 


n.a. 






Marsh Island 


USA 


M 


23 


13 


3 


10 


13 


0.17 




Pyrenees 


France 


M 


19 


13 


5 


6 


13 


0.92 




Redwater 


Canada 


M 


12 


11 


3 


8 


10 


-0.06 




Baikalsee 


Russia 


M 


16 


13 


10 


3 


12 


0.27 




Bialowieza 


Poland 


M 


60 


40 


10 


30 


40 


2.12*** 




Noonday 


USA 


M 


22 


15 


8 


7 


14 


0.11 




Grass Creek 


USA 


M 


28 


22 


11 


11 


21 


-0.22 




Rothwald 


Austria 


M 


43 


23 


12 


11 


23 


0.70* 




Derberence 


Switzerland 


M 


37 


26 


11 


15 


26 


1 .66*** 




Valcartier 


Canada 


M 


35 


30 


17 


13 


28 


0.48 




Bodmeren 


Switzerland 


M 


90 


54 


23 


31 


44 


0.58 
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Table 2 Mating type ratios for populations of species belonging to the Phialocephala fortinii s.l. - Acephala applanata 
species complex (PAC) (Continued) 





Scatle 


Switzerland 


M 


118 


74 


42 


32 




71 


0.71** 




Alpes Maritimes 


France 


M 


125 


75 


46 


29 




75 


0.58* 




Campolino 


Italia 


M 


138 


79 


43 


36 




86 


0.41 




Kevo 


Finland 


M 


138 


87 


44 


43 




72 


0.52* 


P. fortinii 
s.s. 


Glentanar 


Scottland 


M 


10 


5 


3 


2 










Bodmeren 


Switzerland 


M 


10 


8 


4 


4 










Marsh Island 


USA 


M 


17 


10 


5 


5 




10 


-0.17 




Noonday 


USA 


M 


11 


7 


2 


5 










Hill of Fare 


Scottland 


M 


23 


15 


9 


6 




13 


-0.07 




Valcartier 


Canada 


M 


23 


17 


8 


9 




15 


0.37 




Etang de la 


Switzerland 


M 


24 


18 


10 


8 




23 


0.29 


























DldlUWIcZd 


Pnbnrl 
rUldl IU 


M 


Q9 


j i 


9^; 


1 n 




j i 


0 f)^ 




La Facie 


France 


M 


1 09 
I uz 


97 
j / 


1 n 
I O 


21 




99 

JO 


0 99 

U.Z3 




Le Pirou 


France 


M 


1 3 


1 1 


4 


7 




1 1 


0.22 




I oLI IUI 1 IUI IUI d 


Ul\l dll lc 


M 


76 


41 


28 


1 3 


-X- 


41 


0 25 






Pi n 1 a nH 
ril lldl IU 


M 


1 00 


54 


24 


30 




54 


0 28 


LjrO 


f9lnhal 


C-,\nhz\ 
\3\Uua\ 


hA 


ZL9 


90 
zu 


Q 


1 1 










olUUdl 


olUUdl 


M 


1 n 




4 


1 


n.a. 






Ljr 1 U 


olUUdl 


olUUdl 


M 


DU 




4 


1 


n.a. 






CSP1 1 






M 


20 


1 1 


7 


4 




1 1 


0.59 




Kevo 


Pi n 1 a nH 
ril lldl IU 


M 


1 2 


7 


4 


3 








Ljr 1 Z 


C roi iv Hi i 

L.I cUX UU Vdl 1 


OWI IZcl Idl IU 


hA 

IVI 


1 ^ 
i j 


7 
/ 


c: 
D 


9 
z 


n.a. 








Chychkan 


Kirgistan 


M 


9 


7 


6 


1 


n.a. 






CSP13 


Creux du Van 


Switzerland 


M 


43 


23 


9 


12 




23 


0.66 




Alaska 


USA 


M 


14 


10 


7 


3 




11 


0.41 


CSP14 


Kevo 


Finland 


M 


16 


14 


4 


10 




14 


-0.18 




Creux du Van 


Switzerland 


M 


21 


8 


6 


2 


n.a. 








ZinaLBaiting 


Switzerland 


M 


72 


26 


14 


12 




25 


-0.19 


CSP15 


Noonday 


USA 


M 


28 


14 


6 


8 




15 


-0.7 


CSP16 


Creux du Van 


Switzerland 


M 


9 


3 


2 


1 


n.a. 






CSP18 


Global 


Global 


M 


20 


4 


1 


3 


n.a. 






CSP19 


Global 


Global 


M 


30 


5 


4 


1 


n.a. 







1 Molecular markers used for clone-correction were 11 single-copy Restriction Fragment Length Polymorphisms (RFLP) or 12 microsatellite loci (M). 

2 n.a. = not available; significance test was not performed due to small sample size; * = significant deviation from a 1:1 ratio. 

3 MLH = Multilocus haplotypes inferred by RFLP or M; MLHs with missing data were excluded to calculate the index of association (l A ). 

4 * = p < 0.05; ** = p < 0.01; *** = p < 0.001 



Analysis of selective pressures based on sequence 
information 

Overall values of oo ranged between 0.26 and 0.46 for 
the three MAT genes and proportions of amino acid 
sites under purifying selection (<d < 1) were larger than 
proportions of sites reported under neutral selection (co 
= 1) or positive, diversifying selection (oo > 1). For all 
three MAT genes, models that allow codons to evolve 
under positive selection (M2a and M8) did not fit the 
data significantly better than models that do not permit 
positive selection (Mia and M7) (Table 4). Although for 
the genes MAT 1-1-3 and MAT 1-2-1 the log likelihood 
values of Model M2a and M8 were higher than the 



associated models Mia and M7, the ^-values of the LRT 
statistic were > 0.05. Nevertheless, for MAT1-1-3 and 
MAT 1-2-1 the LRT test of M0 versus M3 was signifi- 
cant, suggesting variable selection pressure among sites. 
For these two genes a proportion of sites with co > 1 
was found in models M2a and M8, i.e. 7.6% of the sites 
with co = 2.72 for MAT1-1-3 and 3.1% of the sites with 
co = 5.76 for MAT 1-2-1. Up to seven amino acid sites 
were identified as sites of positive selection using the 
BEB analysis in MAT 1-1-3 but only one site had strong 
support with BEB posterior probability > 95% in model 
M8 (see Additional file 3). In MAT1-2-1 five amino acid 
sites were reported under positive selection but had 
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Table 3 Spatial distribution of mating types for selected PAC populations 

Species Study site Number Fraction of grid points Number of grid Number of grid Association chi- significance 1 

of strains with only one strain points with > 1 points and test 

strain identified mating 

types 













MAT1- 
1 


MAT1- 
2 


both 


[Q] 






P. 


Bodmeren 


28 


0.65 


6 


7 


8 


2 


0.33 


18.75 




europoeo 
























Creux du 
Van 


135 


0.26 


35 


21 


15 


1 1 


-0.03 


4.41 






Zurichberg 


129 


0.24 


42 


20 


21 


14 


-0.22 


3.60 


n.s. 




Scatle 


35 


0.21 


11 


2 


6 


6 


0.94 


14.54 






Derberence 


186 


0.19 


51 


9 


30 


24 


-0.01 


5.11 




P. 

helvetica 


Zurichberg 


40 


0.64 


10 


20 


6 


2 


-0.13 


12.16 




P. 

subolpino 


Bodmeren 


80 


0.46 


22 


5 


25 


11 


0.49 


4.92 






Derberence 


36 


0.70 


7 


9 


12 


2 


-0.03 


12.80 






Scatle 


115 


0.30 


30 


17 


11 


15 


0.43 


3.61 


n.s. 


CSP13 


Creux du 
Van 


43 


0.61 


9 


7 


13 


3 


0.25 


12.06 





1 * = significant at p < 0.05; n.s. = not significant 



weak support with BEB posterior probabilities. In 
MAT1-1-1, no proportions of sites with oo > 1 were 
reported. 

Discussion 

Members of the Phialocephala fortinii s.l. - Acephala 
applanata species complex (PAC) have no known sex- 
ual state. In this study, however, footprints of sexual 
reproduction were found in populations from large glo- 
bal samples, providing evidence that sexual reproduction 
may occur in these fungal species. 

MAT locus structure and evolution of homo- and 
heterothallism in the PAC species 

Homothallic fungi can fertilize themselves whereas het- 
erothallic fungi depend on another compatible indivi- 
dual for sexual reproduction to occur. In filamentous 
ascomycetes the MAT locus is the key determinant of 
breeding system and it has been shown that conversions 
between heterothallism and homothallism can be 
achieved by manipulating the MAT locus [41,49]. For 
example, the heterothallic Neurospora crassa was cap- 
able of self-fertilization after a strain was made carrying 
both mating types [50,51]. In contrast, a strain of the 
homothallic Giberella zeae was converted to self-sterility 
after those MA T genes were deleted that are present on 
opposite mating types in other closely related hetero- 
thallic species [49,52]. 

The structure of the MAT locus was mapped for 11 
additional PAC species for which the MAT locus had 
not yet been characterized. All of these species 



possessed either the MAT 1-1 or MAT 1-2 idiomorph 
consistent with a heterothallic organization structure. 
The identified MAT genes MAT 1-1-1 and MAT 1-3-1 or 
MAT 1-2-1 in the respective idiomorphs were of consis- 
tent lengths and their arrangement within the MAT 
locus and orientation with respect to each other was the 
same as found for other PAC species [44] and the other 
closely related helotialean species Pyrenopeziza brassi- 
cae, Oculimacula yallundae and Rhynchosporium secalis 
[45,46,53]. The only PAC species containing a homo- 
thallic MAT locus structure remains A. applanata [44]. 
A controversial topic in fungal biology is whether het- 
erothallism or homothallism represents the ancestral 
state [37,41,49,54,55]. The question has been addressed 
in particular for the ascomycete genera Cochliobolus 
[40], and references therein] and Aspergillus [56], and 
references therein] that comprise heterothallic and 
homothallic species. In the heterothallic Cochliobolus 
species the MA Tlocus structure is conserved as strains 
are either MAT 1-1 or MAT 1-2 and contain the genes 
MAT 1-1-1 or MAT 1-2-1 with the same gene orientation 
within idiomorphs. In contrast, in the homothallic spe- 
cies the MAT locus structure is unique because single 
individuals possess both MAT genes that are either 
fused into a single ORF, are closely linked or likely pre- 
sent on different chromosomes. Due to this variation in 
MA T locus structure and phylogenetic evidence of inde- 
pendent evolution of self-fertility, the likely derived state 
in this genus was proposed to be homothallism [49]. 
Recombination is normally suppressed in the MAT 
locus because the sequences of the two idiomorphs are 
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Table 4 Likelihood ratio tests comparing models of molecular evolution of MAT genes 



Gene 


Analysis 




df 1 


2AV 


p- 
value 


d N /d s 


Parameters 




MAT1-1- 
7 


MO: one ratio 


-2205.4733 


4 


5.9586 


0.2023 


0.4280 


co = 0.4280 






M3: discrete (« = 3) 


-2202.4940 








0.4462 


Po = 0.3732, co 0 
0.7052 


= 0.7052; Pt = 0.3932, <n, = 0.1 117; p 2 = 0.2336, co 2 = 




M1a: nearly neutral 


-2202.6642 


2 


0.0000 


1 .0000 


0.4599 


Po = 0.6556, co 0 


= 0.1761; Pt = 0.3440, = 1 




M2a: positive 
selection 


-2202.6642 








0.4599 


Po = 0.6556, co 0 


= 0.1761; Pt = 0.2332, ©, = 1; p 2 = 0.1109, co 2 = 1 




M7: p 


-2202.5615 


2 


0.0000 


1 .0000 


0.4507 


Po = 0.4117; q = 


= 0.5019 




M8: (3 & co > 1 


-2202.5615 








0.4507 


Po = 0.1000 p = 


0.4117 q = 0.5019; {p ] = 0.00001) co = 1.0000 


MAT1-1- 

0 

J 


MO: one ratio 


-1706.3379 


4 


20.5175 


0.0004* 


0.295 


co = 0.2951 




M3: discrete (k = 3) 


-1696.0792 








0.3103 


Po = 0.0007, co 0 
2.7212 


= 0.1 116; p! = 0.9232, <o, = 0.1 1 1 7; p 2 = 0.0761 , co 2 = 




M1a: nearly neutral 


-1697.3663 


2 


2.5742 


0.2761 


0.2598 


Po = 0.7402, co 0 


= 0; Pt = 0.2598, = 1 




M2a: positive 
selection 


-1696.0792 








0.3103 


Po = 0.9239, co 0 


= 0.1117; Pt = 0, cot = 1; p 2 = 0.0761, co 2 = 2.7212 




M7: p 


-1697.6108 


2 


3.0594 


0.2166 


0.3000 


p = 0.0050; q = 


0.0119 




M8: p & co > 1 


-1696.0811 








0.3104 


Po = 0.9244 p = 


12.5791 q = 99.0000; fa = 0.0756) co = 2.7290 


MAT1-2- 
1 


MO: one ratio 


-2133.8605 


4 


11.1677 


0.0247* 


0.3714 


co = 0.37135 






M3: discrete (« = 3) 


-2128.2767 








0.4212 


Po = 0.3531, co 0 
5.7427 


= 0.2495; = 0.6157, <n, = 0.2496; p 2 = 0.0312, co 2 = 




M1a: nearly neutral 


-2129.6740 


2 


2.7948 


0.2472 


0.3528 


Po = 0.6472, co 0 


= 0; Pt = 0.3528, = 1 




M2a: positive 
selection 


-2128.2767 








0.4212 


Po = 0.9688, coo 


= 0.2496; p ] = 0, ©! = 1; p 2 = 0.0312, co 2 = 5.7428 




M7: p 


-2129.8169 


2 


3.0800 


0.2144 


0.3387 


p = 0.0170; q = 


0.0321 




M8: p & co > 1 


-2128.2769 








0.4212 


Po = 0.9690 p = 


33.0717 q = 99.0000; fa = 0.0310) co = 5.7612 



Models of neutral evolution (M1a and M7) are tested against models that incorporate positive selection (M2a and M8). Comparison of MO vs. M3 tests variable 
dN/dS rate ratios among sites. Model M3 classifies sites in the sequence into discrete classes (« = 3) with both the dN/dS rate ratio ro 0 , © q , co 3 and the 
proportions p 0 , p 2 and p 3 estimated from the data. For models M1a and M2a, values of p 0 , p 1f and p 2 , are the proportions of sites inferred to be evolving under 
purifying selection, neutral evolution, and positive selection, respectively, and co 0 , ©1 and co 2 are their corresponding dN/dS rate ratios. For models M7 and M8, 
the p distribution, p (p, q), describes the distribution of the dN/dS rate ratio between zero and one, and p 0 is the proportion of sites within this distribution. In 
these models, p q is the proportion of sites inferred to be under positive selection, and co is the dN/dS rate ratio for those sites. 



significantly diverged [57,58]. However, small identity 
island of 8 and 9 nucleotides were identified within 
MA T genes of the heterothallic Cochliobolus heterostro- 
phus and rare recombination events between these small 
identity islands were suggested to be the likely mechan- 
ism for the conversion of reproductive modes in ances- 
tral lineages [49]. In contrast, the prevalence of 
homothallic species, phylogenetic analyses and compari- 
sons of genome sequences suggest that heterothallism is 
the derived state in the genus Aspergillus [59] albeit 
additional characterizations of MAT loci of supposedly 
asexual Aspergillus species may alter this view [56]. 

For PAC species both evolutionary trajectories of 
reproductive modes from homothallism to heterothal- 
lism and vice versa are conceivable. Because MAT genes 
control reproductive processes the information gained 
by comparing their DNA sequences might allow recon- 
structing the evolutionary history of a species complex 
and indicate mechanisms that led to changes in repro- 
ductive modes [49]. The basal position of A. applanata 



in the MAT gene phylogenies suggests that the MAT 
genes of A. applanata are older than the MA T genes of 
the remaining PAC species which are grouped in clus- 
ters, and thus are more recently derived and share a 
common ancestor (see Additional file 4). In fact, this 
finding is congruent with population genetic data and 
other phylogenetic data based on DNA sequences of 
housekeeping genes and RFLP loci that showed that A. 
applanata is the most diverging species in the PAC 
[60]. It is possible that a switch in reproductive mode of 
a homothallic ancestor was the source of the strong 
divergence between A. applanata and other PAC species 
and that thus, A. applanata represents an ancient line- 
age within the PAC. In populations of an ancestral 
homothallic species a structural change in the MAT 
locus, e.g. after a deletion of a MAT gene could have 
provoked a switch to heterothallism and created barriers 
to self-mating. These genetic barriers could have 
increased and maintained genetic diversity, eventually 
leading to speciation and radiation that gave rise to the 
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current heterothallic species. The likely mechanisms that 
have initiated such a switch in the PAC might be 
explained by re-arrangements among idiomorphs and 
the influence of transposable elements on recombination 
events in the MAT locus of A. applanata [44]. On the 
other hand the heterothallic PAC species possess a con- 
served MA T locus structure with congruent MA T gene 
orientation and position within the MAT locus. This 
organization is also consistent with the MAT locus 
structure found in closely related helotialean species 
that are heterothallic [45,46,53]. Thus, the prevalence of 
structurally heterothallic species in the PAC and their 
conserved MAT locus organization are arguments in 
favor of heterothallism as ancestral breeding system. 

Mating type distribution and gametic equilibrium 

If sexual reproduction plays an important role in the 
life cycle of a fungal species, natural populations would 
be expected to show following characteristics [16,19]: 
(i) 1:1 ratio of mating types because frequency-depen- 
dent selection acts on the rare mating type that has 
the higher chance for mating, and (ii) gametic equili- 
brium at unlinked loci because under random mating 
sexual recombination will reduce non-random associa- 
tion of alleles. Both criteria were met for the majority 
of PAC populations. The 6 out of 52 populations that 
displayed gametic disequilibrium were restricted to the 
species P. subalpina. Independently of the reproductive 
mode several factors can cause non-random associa- 
tion of alleles that result in gametic disequilibrium. 
These factors include high frequency of rare alleles, 
genetic drift, population admixture and selection 
[16,61]. We therefore re-analyzed these collections of 
P. subalpina after pooling rare alleles but the results 
did not change (data not shown). We then tested 
whether traces of population admixture could be iden- 
tified in the two populations using the program 
STRUCTURE [62] and evidence for admixture was 
found in two populations. Moreover, the strains col- 
lected from Bialowieza could be allocated with high 
posterior probabilities to two subpopulations and the 
index of association I A did not significantly deviate 
from zero in both of these subpopulations (data not 
shown). 

The spatial distribution of mating type within study 
sites is often neglected in studies aiming at assessing the 
reproductive biology of fungal species. In this study 
opposite mating types were often collected at the same 
grid points within study sites. Moreover, the measures 
of spatial association Q were predominantly neutral to 
strongly positive, meaning that for several populations 
the occurrence of both mating types at the same grid 
point was higher than expected by chance alone. Thus, 
the physical proximity of compatible mating types in the 



study sites meets another prerequisite for successful sex- 
ual reproduction in the PAC. 

Selection acting on mating type genes in PAC 

If MAT genes play an important role in the initiation 
and further steps of sexual reproduction [39,40] they are 
expected to be functionally conserved in sexual species. 
All PAC species possessed structurally identical idio- 
morphs and translation of the MAT genes suggested 
that they encode functional proteins. In A. applanata 
the insertion of a transposable element in the 3' end of 
MAT 1-2-1 caused a truncation of 30 amino acids in the 
translated protein [44]. Moreover, insertions of 3 and 9 
bp were found in MAT1-1-1 of A. applanata and in 
MAT 1-2 deletions as long as 18 bp were found for cryp- 
tic species CSP12, CSP14, P. letzii, P< europaea and A. 
applanata. However these indels did not result in stop 
codons or frame shifts that could affect functionality. 

The molecular analyses revealed that all three MAT 
genes appear to be under strong purifying selection 
because sites under purifying selection were predomi- 
nant. The overall values of <d obtained ranged between 
0.26-0.46, and the proportions of sites under purifying 
selection (co < 1) were much larger than the fraction 
evolving neutrally or under positive selection accounting 
for at least 65% of the sites. In only few studies the 
selective pressures on MAT genes were analyzed and 
similar overall values of co were found ranging between 
0.14 and 0.49 in sexual or presumed sexual fungi 
[22,23,25,26]. The predicted effect of long-term asexual - 
ity is the decay of genes specific for sex and recombina- 
tion due to the relaxation of selective constraints [63] 
but the results of our study suggest that MA T genes in 
PAC species are highly conserved and functional. 

Rapid divergence of reproductive proteins, due to 
positive selection and adaptive evolution, is believed to 
be important in the formation of reproductive barriers 
leading to speciation [64-66]. Competition for gamete 
recognition between individuals might explain the rapid 
evolution of reproductive proteins [67]. High levels of 
interspecific polymorphism have been found in MAT 
genes [22,68]. In our study the amount of polymorphism 
between the structurally heterothallic PAC species was 
also high and ranged between 8.7 and 9.4%. This was in 
the range found for non-coding RFLP loci and approxi- 
mately 2.5-fold more than in levels found in the house- 
keeping gene beta-tubulin [60]. However, statistical 
evidence of positive selection of MAT gene, i.e. high 
levels of non-synonymous substitutions leading to 
changes in MAT proteins, has only been shown in the 
genus Neurospora [25]. Our results showed that models 
of evolution that allow positive selection (co > 1) were 
not more likely than those incorporating neutral evolu- 
tion. Nevertheless, the comparison of model M0 vs. M3 



Zaffarano et al. BMC Evolutionary Biology 201 1, 1 1:282 
http://www.biomedcentral.eom/1 471 -21 48/1 1 /282 



Page 10 of 14 



tested significant for variable pressure among sites in 
MAT 1-1-3 and MAT 1-2-1 and up to seven sites evolving 
under positive selection were identified in these genes, 
although only one was significant in the BEB analysis. It 
might be assumed that these and many more sites 
evolve under positive selection in MAT genes but that, 
in a strong background of purifying selection, they 
remain uncovered by the currently available analysis 
methods. However, whether the putative protein altera- 
tions at such positively selected sites have a biological 
significance is unknown. 

Have PAC species a sexual life cycle? 

In the present study we combined different indirect 
approaches to analyze the importance of sexual repro- 
duction in the PAC. The population analyses were in 
accordance with random mating and sexual recombina- 
tion in abundant collections from different ecosystems. 
Moreover, the presence of a conserved MAT locus in all 
PAC species suggests that MAT genes are functional 
and are still involved in mating processes. Yet, although 
the data does not allow rejecting the hypothesis of sex- 
ual recombination, the confirmation of the teleomorph 
in vitro is pending. Our first attempts to induce apothe- 
cia formation of PAC species failed when we tried to 
cross opposite mating types under laboratory conditions 
[44]. However, in vitro induction of the teleomorph is 
often difficult and time consuming and some persever- 
ance is required [17,18]. Nevertheless, we successfully 
induced teleomorph formation of Phaeomollisia piceae, 
the closest known relative of the PAC [69]. P. piceae 
forms small darkly pigmented apothecia of less than 800 
[im diameter. Considering that PAC species might form 
similar apothecia in the soil or on roots it seems likely 
that the teleomorph has been overlooked in field 
studies. 

Conclusions 

In this study we evaluated the importance of sexual 
reproduction in the PAC, the dominant dark septate 
root endophytes of woody plant species, and showed 
that the signature of sex is present in all of these puta- 
tively asexual species. The MAT locus resembles those 
of closely related heterothallic ascomycetes, is con- 
served and purifying selection likely preserves func- 
tionality of its genes. The hypothesis of sex occurring 
in the PAC cannot be rejected based on the random 
association of alleles at multiple loci, the equal mating 
type frequencies and the spatial distribution of oppo- 
site mating types in populations from different ecologi- 
cal habitats, hosts and continents. We believe that 
further field studies and in vitro crosses will finally 
lead to the discovery of the sexual state, which is the 
missing link to better understand how reproductive 



biology shaped the evolutionary history and influences 
the ecological role of the PAC. 

Methods 

Collection of PAC species and mating type screening 

A total of 3, 639 strains derived from 33 study sites 
were included in the mating type screening (Table 2). 
These strains represent 19 PAC [32] species eight of 
which are formally described [28,29]. Several classes of 
molecular markers were developed for PAC species 
assignment including PCR fingerprints, single-copy 
restriction fragment length polymorphisms (RFLP), 
multi-locus DNA sequences, and microsatellites. Each of 
these molecular markers supported the delineation of 
multiple species in this complex, with concordant cryp- 
tic species defined by all markers [60,70]. PAC species 
communities were either collected on grid nets (196 m 2 , 
74 grid points) [33] or on transects (50 m, 11 grid 
points) [32] from Europe (number of study sites n = 
19), North America (n = 10) and Asia (n = 4). Collec- 
tions included both undisturbed and managed forests. 

Multiplex PCR 

A multiplex PCR method was developed to discriminate 
between isolates carrying either the MAT 1-1 or the 
MAT 1-2 idiomorph based on available sequence data of 
the MA T locus of seven structurally heterothallic PAC 
species [44]. The primer Pf_HMG_R.03 (5'- 
TCCTCGAACCGGTGCCACGAATGACTCCGA-3') 
was designed to match a flank of the idiomorphic region 
that is conserved between the two mating types of these 
species. The MAT 1-1 specific primer (Pf_HMG_F4: 5'- 
AGCTGAGCTCGAGACCTCGTTCGC-3') and a 
MAT 1-2 specific primer (Pf_MATl-lFlc: 5'-CTTGTA 
CGGTCCGGCAATCCACA-3') were designed to gener- 
ate DNA fragments of distinct length in combination 
with Pf_HMG_R.03 that could be easily differentiated 
on agarose gels. All PCRs were performed on a Biome- 
tra Tl thermal cycler in a 10 \A reaction volume con- 
taining approximately 2 ng template DNA, 5 pmol of 
each primer, 50 mM KC1, 10 mM Tris-HCl, 1, 5 mM 
MgCl2, 200 uM dNTPs (Promega AG, Dubendorf) and 
0.2 U Taq polymerase (Promega AG, Dubendorf). 
Cycling conditions during PCR were 2 min at 94°C fol- 
lowed by 29 cycles of denaturation for 20 s at 94°C, 
annealing for 30 s at 60°C and extension for 60 s at 72° 
C, followed by a final extension step of 5 min at 72°C. 

Mating type frequencies in PAC populations 

The number of strains from a single study site and 
belonging to the same PAC species is here defined as a 
population. Before assessing significance of deviation 
from a 1:1 mating type ratio in a population, each popu- 
lation was clone-corrected using 12 microsatellite loci 
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[70] or 11 single-copy RFLP markers [33], meaning that 
only one representative of each multi-locus genotype 
was included. Clone-correction is necessary because the 
same fungal thallus is found on several grid points. Five 
species [cryptic species (CSP) 8, 9, 10, 18 and 19] 
occurred at exceptionally low frequencies in the collec- 
tion sites and represented less than 3% of all collected 
strains. For these species, the distribution of mating 
types was assessed after pooling all strains of a species, 
irrespective of the study site from which they were iso- 
lated. Significance was tested using the exact binomial 
test. In addition, Bonferroni corrections were performed 
for multiple comparisons. 

Because species abundance distributions within local 
PAC communities follow a hyperbolic distribution with 
a few very abundant species and many "rare" species 
[48], clone-corrected datasets for such rare species 
resulted in a low number of distinct genotypes (< 10). 
In these cases, only the frequency of mating types was 
recorded without testing for significance. 

Analysis of gametic disequilibrium 

Significant deviations from multi-locus gametic equili- 
brium in a population were tested using the index of 
association I A and its variance according to Maynard 
Smith et al. [47]. All calculations of I A were based on 
clone-corrected datasets and multi-locus genotypes with 
missing data were excluded from the analysis. 

Spatial distribution of mating types in selected study 
sites 

Spatial distribution of mating types was analyzed for five 
Swiss study sites. Only the two most abundant PAC spe- 
cies per community were included in this analysis. For 
each of these species, the number of grid points was 
recorded where either MAT 1-1 or MAT 1-2 or both 
mating types were found. A contingency table test for 
independence of the spatial distribution of the mating 
types was performed and a measure of association 
between MAT 1-1 and MAT 1-2 (Q) was calculated as 
described in Pielou [71]. 

Studying the presence of both mating types in single 
isolates 

In few cases, no amplification product was obtained in 
some strains, or both the MAT1-1 and MAT1-2 speci- 
fic fragments were amplified in single strains. In order 
to study the reason for these observations, ten new 
single-hyphal tip cultures were prepared from the ori- 
ginal slants of the three strains having displayed both 
PCR products and of the five strains from which no 
PCR product was obtained. After harvesting and 
extracting DNA of mycelia from both the original 
slants and the newly prepared single-hyphal tip 



cultures, the mating types were re-determined using 
the multiplex PCR. 

Sequencing the MAT locus 

The MAT locus was sequenced for 11 PAC species and 
compared to the MAT locus structure of 8 previously 
sequenced PAC species [44]. For each of these species 
the MAT locus of both the idiomorph MAT 1-1 and 
MAT 1-2 was sequenced (Table 1). The detailed sequen- 
cing strategy and the sequencing primers used are given 
in the Additional files 5 &6, and the MAT locus was 
annotated as described in [44]. PCR fragments were 
directly purified using an "ExoSap" protocol [70]. Cycle 
sequencing was performed with the Big-Dye v3.1 kit 
(Applied Biosystems) in a 10 |il reaction volume using 1 
\i\ purified PCR amplification product, 0.5 [i\ BigDye 
v3.1, 1.9 ul 5 x Buffer, 5.6 [i\ ddH20 and 1 [i\ of the 
sequencing primer (10 \iM). Cycle sequencing reactions 
were performed on a Biometra Tl thermal cycler with 
the following running conditions: 60 s at 96°C followed 
by 55 cycles of 10 s at 95°C, 5 s at 50°C and 4 min at 
60°C. Dye-labeled fragments were cleaned using BigDye 
Xterminator Purification kit following manufacturer's 
instructions (Applied Biosystems). Samples were run on 
an ABI 3130x1 DNA Analyser (Applied Biosystems). 

Analysis of selective pressures 

In order to assess what kind of selective pressure is act- 
ing on any of the three MA T genes a maximum likeli- 
hood (ML) method implemented in the program 
CODEML in the PAML 4.1 software package [72] was 
used to compare the rate of non-synonymous substitu- 
tions with the rate of synonymous substitutions (d N /d s 
= <d). Under neutrality <d is expected to be 1. When 
amino acid changes are favored, indicative of positive 
selection, go is expected to be > 1, whereas under purify- 
ing selection amino acid changes are prevented and co is 
expected to be < 1. Positive selection of small number 
of codons in the MA T genes may be masked when the 
entire gene is strongly affected by purifying selection. 
Therefore, a codon by codon approach was used based 
on the ML method implemented in CODEML that 
allows site by site identification of particular codons that 
have been evolving under repeated and strong positive 
selection [72]. The three MAT genes were analyzed 
separately and a gene tree describing the phylogenetic 
relationships of all taxa studied (one strain per species) 
was constructed using the ML method implemented in 
the program treefinder [73] using the best-fitting muta- 
tion model. The inferred trees served as the basis for 
the implementation of the ML methods in the CODML 
(see Additional file 4). 

Three sets of models, commonly applied to test 
hypotheses of selection, were used. Following the 
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suggestions of Yang [72] the site model pairs that 
appear to be particularly useful for real data analysis, are 
the Mia (nearly neutral) versus M2a (selection) and M7 
(neutral) versus M8 (selection). However, model MO 
(one-ratio) was also compared versus M3 (discrete) in 
order to see if the selective pressure is variable among 
sites. Because the comparison of M7 versus M8 is less 
conservative, this set of models may indicate positive 
selection even when none is detected by the Mia vs. 
M2a comparison. 

The strength of positive selection was calculated using 
the Likelihood Ratio Test (LRT) [72] implemented in 
PAML by comparing twice the log likelihood difference 
in a chi-square test. The degree of freedom for the test 
statistic is determined by the difference in estimated 
parameters between the models of selective pressures 
being compared, i.e. four [MO (one-ratio) vs. M3 (dis- 
crete)] or two [Mia (nearly neutral) vs. M2a (selection), 
and M8 (selection) vs. M7 (neutral)] degrees of freedom 
[72]. Codons that are identified as having evolved under 
positive selection have high posterior probabilities (p > 
0.95). Posterior probabilities for such sites were esti- 
mated based on Bayes Empirical Bayes (BEB) analysis 
[74]. 

The MAT gene sequence of the homothallic A. appla- 
nata species was not included in this analysis because it 
would have inflated the alignments with deletions. For 
example, the MAT1-2-1 sequence of A. applanata con- 
tained a large deletion of 90 bp at its 3' end, due to the 
insertion of a transposable element in the MA T locus of 
the species [44]. Gaps in alignments cannot be handled 
by CODEML and the program automatically removes 
such sites. Thus, in order to prevent the loss of informa- 
tive sites in the sequence of the other species, the A. 
applanata sequences were removed from the alignment 
prior to the analyses. 

Additional material 



Additional file 6: Sequencing primers. Sequences of primers used to 
sequence the complete MAT idiomoprh of PAC species. 
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complex. 

Acknowledgements 

The Genetic Diversity Center (GDC, ETH Zurich) provided facilities for 
collecting the sequence data. This work was associated with the GEDIHAP 
project of the Competence Center Environment and Sustainability (CCES) of 
the ETH Domain, Zurich, Switzerland. 

Author details 

institute of Integrative Biology (IBZ), Forest Pathology and Dendrology, ETH 
Zurich, 8092 Zurich, Switzerland. 2 Microsynth AG, Schutzenstrasse 15, 9436 
Balgach, Switzerland. 

Authors' contributions 

PZ participated in the design of the study, conducted the analysis of 
selective pressures, drafted and wrote the manuscript. VQ participated in the 
molecular work and conducted microsatellite analysis. AD carried out 
molecular work and participated in the sequence alignment. CG designed 
and coordinated the study, performed statistical analysis and helped to draft 
the manuscript. All authors read and approved the final manuscript. 

Received: 5 May 201 1 Accepted: 30 September 201 1 
Published: 30 September 2011 

References 

1. Bell G: The Masterpiece of nature: the evolution and genetics of 
sexuality. London Canberra: Croom Helm; 1982. 

2. Kondrashov AS: Classification of hypothesis on the advantage of 
aphimixis. J Hered 1993, 84(5):372-387. 

3. West S, Lively CM, Read AF: A pluralist approach to sex and 
recombination. J Evolution Biol 1999, 1 2(6):1 003-1 01 2. 

4. Agrawal AF: Evolution of sex: Why do organisms shuffle their genotypes? 
Curr Biol 2006, 1 6(1 7):R696-R704. 

5. Otto SP: The Evolutionary Enigma of Sex. Amer Natur 2009, 1 7451-S14. 

6. Meirmans S, Strand R: Why Are There So Many Theories for Sex, and 
What Do We Do with Them? J Hered 2010, 101:S3-S12. 

7. Crow JF, Kimura M: Evolution in sexual and asexual populations. Amer 
Natur 1965, 99(909):439-450. 

8. Smith JM: The evolution of sex Cambridge: Cambridge University Press; 1978. 

9. Charlesworth B: The Cost of Sex in Relation to Mating System. J Theor Biol 
1980, 84(4):655-671. 

10. Hoekstra RF, Vanloo EN: The Cost of Sex in Hermaphrodite Populations 
with Variation in Functional Sex. J Theor Biol 1986, 1 22(4):441 -452. 

11. Kondrashov AS: Deleterious Mutations and the Evolution of Sexual 
Reproduction. Nature 1988, 336(61 98):435-440. 

12. Aanen D, Hoekstra RF: Why sex is good: On fungi and beyond. In Sex in 
Fungi Molecular Determination and Evolutionary Implications. Edited by: 
Heitman J, Kronstad JW, Taylor JW. Casselton LA. Washington, DC: ASM 
Press; 2007:527-534. 

13. Kirk PM, Cannon PF, David JC, Stalpers JA, eds: Dictionary of the fungi. 
Oxon, UK: CAB International;, 9 2001. 

14. Dyer PS, Ingram DS, Johnstone K: The Control of Sexual Morphogenesis in 
the Ascomycotina. Biol Rev 1992, 67(4):421-458. 

15. Kiick U, Poeggeler S: Cryptic sex in fungi. Fung Biol Rev 2009, 23:86-90. 

16. Milgroom MG: Recombination and the multilocus structure of fungal 
populations. Annu Rev Phytopathol 1996, 34:457-477. 

17. Magee PT, Magee BB: Through a glass opaquely: the biological 
significance of mating in Candida albicans. Curr Opin Microbiol 2004, 
7(6):661-665. 

18. O'Gorman CM, Fuller HT, Dyer PS: Discovery of a sexual cycle in the 
opportunistic fungal pathogen Aspergillus fumigatus. Nature 2009, 

457(7228):471-U475. 



Additional file 1: Testing the presence of both idiomorphs in 
selected PAC strains. Testing the presence of both idiomorphs, or the 
absence of MAT genes in selected PAC strains. 

Additional file 2: Spatial distribution of mating types. Example of the 
spatial distribution of mating types for a population of Phialocephala 
subalpina and Phialocephala turicensis. 

Additional file 3: Sites in MAT genes identified to evolve under 
positive selection. Sites under positive selection were identified using 
the Bayes Empirical Bayes analysis under different CODML site models in 
the PAML software package. 

Additional file 4: Phylogenies inferred from MAT genes. Phylogenies 
inferred for the three genes MAT1-1-1, MAT! -1-3, and MAT1-2-I. The 
sequence of Acephala applanata served as outgroup in all three analysis. 

Additional file 5: Sequencing strategy for MAT idiomorphs. 

Sequencing strategy applied to sequence the complete MAT idiomorphs 
for PAC species. Sequences of primers are given in the Additional file 5. 



Zaffarano et al. BMC Evolutionary Biology 201 1, 1 1:282 
http://www.biomedcentral.eom/1 471 -21 48/1 1 /282 



Page 13 of 14 



19. Xu J: Fundamentals of fungal molecular population genetic analyses. 

Curr Issues Mol Biol 2006, 8:75-89. 

20. Taylor JW, Jacobson DJ, Fisher MC: The evolution of asexual fungi: 
reproduction, speciation and classification. Annu Rev Phytopathol 1999, 
37:197-246. 

21. Paoletti M, Rydholm C, Schwier EU, Anderson MJ, Szakacs G, Lutzoni F, 
Debeaupuis JP, Latge JP, Denning DW, Dyer PS: Evidence for sexuality in 
the opportunistic fungal pathogen Aspergillus fumigatus. Curr Biol 2005, 
15(13):! 242-1 248. 

22. O'Donnell K, Ward TJ, Geiser DM, Kistler HC, Aoki T: Genealogical 
concordance between the mating type locus and seven other nuclear 
genes supports formal recognition of nine phylogenetically distinct 
species within the Fusarium grominearum clade. Fung Genet Biol 2004, 
41(6):600-623. 

23. Devier B, Aguileta G, Hood ME, Giraud T: Ancient Trans-specific 
Polymorphism at Pheromone Receptor Genes in Basidiomycetes. 
Genetics 2009, 181(1):209-223. 

24. Karlsson M, Nygren K, Johannesson H: The evolution of the pheromonal 
signal system and its potential role for reproductive isolation in 
heterothallic Neurospora. Mol Biol Evol 2008, 25(1 ):1 68-1 78. 

25. Wik L, Karlsson M, Johannesson H: The evolutionary trajectory of the 
mating-type (mat) genes in Neurospora relates to reproductive behavior 
of taxa. BMC Evol Biol 2008, 8. 

26. Lopez-Villavicencio M, Aguileta G, Giraud T, de Vienne DM, Lacoste S, 
Couloux A, Dupont J: Sex in Penicillium: Combined phylogenetic and 
experimental approaches. Fung Genet Biol 2010, 47(8):693-706. 

27. Grunig CR, Queloz V, Sieber TN, Holdenrieder O: Dark septate endophytes 
(DSE) of the Phiolocephola fortinii s.l. - Acepholo opplonato species 
complex in tree roots - classification, population biology and ecology. 
Botany 2008, 86(1 2):1 355-1 369. 

28. Grunig CR, Duo A, Sieber TN, Holdenrieder O: Assignment of species rank 
to six reproductively isolated cryptic species of the Phialocephala fortinii 
s.\.-Acephala opplanata species complex. Mycologio 2008, 100(1):47-67. 

29. Grunig CR, Sieber TN: Molecular and phenotypic description of the 
widespread root symbiont Acepholo opplonoto gen. et sp. nov., formerly 
known as dark septate endophyte Type 1 . Mycologio 2005, 97(3):628-640. 

30. Zhang C, Yin L, Dai S: Diversity of root-associated fungal endophytes in 
Rhododendron fortunei in subtropical forests of China. Mycorrhiza 2009, 
19(6):4 17-423. 

31. Piercey MM, Graham SW, Currah RS: Patterns of genetic variation in 
Phiolocepholo fortinii across a broad latitudinal transect in Canada. Mycol 
Res 2004, 108(8):955-964. 

32. Queloz V, Sieber TN, Holdenrieder , McDonald BA, Grunig CR: No 
biogeographical pattern for a root-associated fungal species complex. 
Global Ecol Biogeog 201 1, 20(1 ):1 60-1 69. 

33. Grunig CR, McDonald BA, Sieber TN, Rogers SO, Holdenrieder O: Evidence 
for subdivision of the root-endophyte Phiolocepholo fortinii into cryptic 
species and recombination within species. Fung Genet Biol 2004, 
41(7):676-687. 

34. McGill BJ, Etienne RS, Gray JS, Alonso D, Anderson MJ, Benecha HK, 
Dornelas M, Enquist BJ, Green JL, He FL, et al: Species abundance 
distributions: moving beyond single prediction theories to integration 
within an ecological framework. Ecol Letters 2007, 1 0(1 0):995-1 01 5. 

35. de Wit R, Bouvier T: 'Everything is everywhere, but, the environment 
selects'; what did Baas Becking and Beijerinck really say? Environ 
Microbiol 2006, 8(4):755-758. 

36. Heitman J, Kronstad JW, Taylor JW, Casselton LA, eds: Sex in Fungi. 
Molecular Determination and Evolutionary Implications. Washington, D.C.: 
ASM Press; 2007. 

37. Coppin E, Debuchy R, Arnaise S, Picard M: Mating types and sexual 
development in filamentous ascomycetes. Microbiol Mol Biol Rev 1997, 
61 (4):41 1-428. 

38. Casselton LA: Mate recognition in fungi. Heredity 2002, 88:142-147. 

39. Stanton BC, Hull CM: Mating-type locus control of cell identity. In Sex in 

Fungi Molecular Determination and Evolutionary Implications. Edited by: 
Heitman J, Kronstad JW, Taylor JW. Casselton LA. Washington, DC: ASM 
Press; 2007:59-73. 

40. Turgeon BG, Debuchy R: Cochliobolus and Podosporo: mechanisms of sex 
determination and the evolution of reproductive lifestyle. In Sex in Fungi 
Molecular Determination and Evolutionary Implications. Edited by: Heitman J, 



Kronstad JW, Taylor JW. Casselton LA. Washington, DC: ASM Press; 
2007:93-121. 

41. Lin X, Heitman J: Mechanisms of homothallism in fungi and transitions 
between heterothallism and homothallism. In Sex in Fungi Molecular 
Determination and Evolutionary Implications. Edited by: Heitman J, Kronstad 
JW, Taylor JW. Casselton LA. Washington, DC: ASM Press; 2007:35-57. 

42. Debuchy R, Turgeon BG: Mating-type structure, evolution, and function in 
Euascomycetes. In The Mycota: Growth, differentiation and sexuality. Volume 
I. Edited by: Kues U, Fischer R. Berlin Springer; 2006:293-323. 

43. Turgeon G, Yoder OC: Proposed nomenclature for mating type genes of 
filamentous ascomycetes. Fung Genet Biol 2000, 31 (1 ):1 -5. 

44. Zaffarano PL, Duo A, Grunig CR: Characterization of the mating type 
{MAT) locus in the Phiolocepholo fortinii s.l. -Acepholo opplonoto species 
complex. Fungal Genet Biol 2010, 47(9):761-772. 

45. Singh G, Ashby Alison M: Cloning of the mating type loci from 
Pyrenopezizo brossicoe reveals the presence of a novel mating type gene 
within a discomycete MAT 1-2 locus encoding a putative 
metallothionein-like protein. Mol Microbiol 1998, 30(4):799-806. 

46. Foster SJ, Fitt BDL: Isolation and characterisation of the mating-type 
{MAT) locus from Rhynchosporium secolis. Curr Genetics 2003, 
44(5):277-286. 

47. Smith JM, Smith NH, Orourke M, Spratt BG: How clonal are bacteria? Proc 
Natl Acad Sci USA 1993, 90(10):4384-4388. 

48. Grunig CR, Duo A, Sieber TN: Population genetic analysis of Phiolocepholo 
fortinii s.l. and Acepholo opplonoto in two undisturbed forests in 
Switzerland and evidence for new cryptic species. Fung Genet Biol 2006, 
43(6)410-421. 

49. Yun S-H, Berbee ML, Yoder OC, Turgeon BG: Evolution of the fungal self- 
fertile reproductive life style from self-sterile ancestors. Proc Natl Acad Sci 
USA 1999, 96(10):5592-5597. 

50. Glass NL, Grotelueschen J, Metzenberg RL: Neurospora crosso A mating- 
type region. Proc Natl Acad Sci USA 1990, 87(1 3):491 2-491 6. 

51. Staben C, Yanofsky C: Neurospora crossa a mating-type region. Proc Natl 
Acad Sci USA 1990, 87(1 3):491 7-4921. 

52. Lee J, Lee T, Lee YW, Yun SH, Turgeon BG: Shifting fungal reproductive 
mode by manipulation of mating type genes: obligatory heterothallism 
of Gibberella zeae. Mol Microbiol 2003, 50(1 ):1 45-1 52. 

53. Singh G, Dyer Paul S, Ashby Alison M: Intra-specific and inter-specific 
conservation of mating-type genes from the discomycete plant- 
pathogenic fungi Pyrenopezizo brassicoe and Tapesio yallundoe. Curr 
Genet 1999, 36(5):290-300. 

54. Nauta MJ, Hoekstra RF: Evolution of reproductive systems in filamentous 
ascomycetes. II. Evolution of hermaphroditism and other reproductive 
strategies. Heredity 1992, 68:537-546. 

55. Butler G: The evolution of MAT: The ascomycetes. In Sex in fungi: 
molecular determination and evolutionary implications. Edited by: Heitman J, 
Kronstad JW, Taylor JW. Casselton LA. Washington, DC: ASM Press; 
2007:3-18. 

56. Dyer PS: Sexual reproduction and significance of MAT in Aspergilli. In Sex 

in Fungi Molecular Determination and Evolutionary Implications. Edited by: 
Heitman J, Kronstad JW, Taylor JW. Casselton LA. Washington, DC: ASM 
Press; 2007:123-142. 

57. Menkis A, Jacobson DJ, Gustafsson T, Johannesson H: The mating-type 
chromosome in the filamentous ascomycete Neurospora tetrosperma 
represents a model for early evolution of sex chromosomes. PLoS 
Genetics 2008, 4(3):e1 000030. 

58. Bergero R, Charlesworth D: The evolution of restricted recombination in 
sex chromosomes. Trends Ecol Evol 2009, 24(2):94-102. 

59. Galagan JE, Calvo SE, Cuomo C, Ma LJ, Wortman JR, Batzoglou S, Lee SI, 
Basturkmen M, Spevak CC, Clutterbuck J, et al: Sequencing of Aspergillus 
nidulons and comparative analysis with A. fumigatus and A. oryzoe. 
Nature 2005, 438(7071 ):1 105-1 1 15. 

60. Grunig CR, Brunner PC, Duo A, Sieber TN: Suitability of methods for 
species recognition in the Phiolocepholo fortinii -Acepholo opplonoto 
species complex using DNA analysis. Fung Genet Biol 2007, 44(8):773-788. 

61. Linde CC, Zala M, McDonald BA: Molecular evidence for recent founder 
populations and human-mediated migration in the barley scald 
pathogen Rhynchosporium secolis. Mol Phylogenet Evol 2009, 51(3):454-464. 

62. Pritchard JK, Stephens M, Donnelly P: Inference of population structure 
using multilocus genotype data. Genetics 2000, 155(2):945-959. 



Zaffarano et al. BMC Evolutionary Biology 201 1, 1 1:282 
http://www.biomedcentral.eom/1 471 -21 48/1 1 /282 



Page 14 of 14 



63. Normark BB, Judson OP, Moran NA: Genomic signatures of ancient 
asexual lineages. Biol J Linnean Soc 2003, 79(1):69-84. 

64. Snook RR, Chapman T, Moore PJ, Wedell N, Crudgington HS: Interactions 
between the sexes: new perspectives on sexual selection and 
reproductive isolation. Evol Ecol 2009, 23(1):71-91. 

65. Turner LM, Hoekstra HE: Causes and consequences of the evolution of 
reproductive proteins. Int J Dev Biol 2008, 52(5-6):769-780. 

66. Clark NL, Aagaard JE, Swanson WJ: Evolution of reproductive proteins 
from animals and plants. Reproduction 2006, 1 31 (1 ):1 1-22. 

67. Palumbi SR: Speciation and the evolution of gamete recognition genes: 
pattern and process. Heredity 2009, 1 02(1 ):66-76. 

68. Poeggeler S: Phylogenetic relationships between mating-type sequences 
from homothallic and heterothallic ascomycetes. Curr Genet 1999, 
36(4)222-231. 

69. Grunig CR, Queloz V, Duo A, Sieber TN: Phylogeny of Phaeomollisia piceae 
gen. sp. nov.: a dark-septate conifer-needle endophyte and its 
relationships to Phialocepholo and Acephala. Mycol Res 2009, 

113(2):207-221. 

70. Queloz V, Duo A, Sieber TN, Grunig CR: Microsatellite size homoplasies 
and null alleles do not affect species diagnosis and population genetic 
analysis in a fungal species complex. Mol Ecol Ressources 2010. 

71. Pielou EC: Association between pairs of species I: Individuals in discrete 
habitable units. Mathematical Ecology New York: John Wiley & Sons; 1977, 
203-239. 

72. Yang ZH: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol 
Evol 2007, 24(8):1 586-1 591. 

73. Jobb G: TREEFINDER version 2010 [http://www.treefinder.de/]. 

74. Yang ZH, Wong WSW, Nielsen R: Bayes empirical Bayes inference of 
amino acid sites under positive selection. Mol Biol Evol 2005, 
22(4):1 107-1 118. 



doi:1 0.1 1 86/1 471 -21 48-1 1 -282 

Cite this article as: Zaffarano et al:. Sex in the PAC: A hidden affair in 
dark septate endophytes? BMC Evolutionary Biology 201 1 1 1:282. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 

Submit your manuscript at (^\ RioMM i r pnt ^\ 

www.biomedcentral.com/submit Blomea ^ enirai 



